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Abstract 

We compute the effect of a non-zero lepton chemical potential on the structure of the three 
flavor Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase of QCD at finite temperature. We show 
that, as in the BCS case, the lepton chemical potential favors two-species color superconductivity 
and disfavors the three species pairing. We stress that this study could be relevant for the cooling 
of a proto-neutron star with a FFLO core, if the temperatures are higher than the un-trapping 
temperature. 
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I. INTRODUCTION 



Investigations on the phase diagram of Quantum-Chromo-Dynamics (QCD) in the high 
quark density and low temperature regime has attracted a lot of interest in the last years; 
apart the purely theoretical speculations, these studies might be relevant for a deeper un- 
derstanding of the physics of compact stellar objects, where cold and dense quark matter 
could be present. In this regime QCD predicts Cooper pairing of quarks, due to the exis- 
tence of an attractive quark interaction in the color antisymmetric channel, see and for 
reviews jsj . Formally, one introduces a bilinear quark expectation value (namely a di-quark 
condensate) in order to describe the collective pairing; since a pair of quarks is not a color 
singlet under SU{3)coior, the condensate spontaneously breaks the color symmetry. This 
phenomenon is similar to the Cooper pairing in ordinary (electromagnetic) BCS supercon- 
ductors where the role of the quarks is played by the electrons and SU{3)coior is replaced 
by the electromagnetic U{1). Because of this analogy the quark condensation phenomenon 
is known as Color Superconductivity. The spectrum of the color superconductive phases is 
usually characterized by gapped fermions, massive gluons and Goldstone bosons related to 
the breaking of some of the global symmetries of the QCD lagrangian. 

At asymptotic densities and zero temperature, the ground state of color superconductivity 
with massless up (n), down {d) and strange (s) quarks is the Color-Flavor-Locking (CFL) 
phase ^. In the intermediate densities regime, as could be found in the interior of a compact 
stellar object, one cannot neglect neither the strange quark mass nor the differences 6fi in 
the quark chemical potentials, induced by /3 equilibrium and electric and color neutrality 
constraints. Therefore in the pre-asymptotic regime the CFL phase could be replaced by 
another ground state. Several ground states have been proposed in the literature as the 
candidates for the "pre-CFL" phases; we recall here the 2SC phase and the gapless 
phases g2SC and gCFL Q (see Q for recent studies). In all the above phases the 
Cooper pairs are characterized by a vanishing total momentum. This allows the whole 
Fermi sphere to partecipate to the pairing. Moreover, the pairs have vanis hing tot al sp in. 
However, the gapless phases present chromo- magnetic instability ^ 11, ll'l. Il3. llJ ] as 
they show imaginary gluon Meissner masses: this should be intimately connected to the 
existence of the gapless modes in these phases. An instability is present also in the 2SC 
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magnetic instability of the 2SC and of the gapless phases). 

For appropriate values of 5yU, it can be advantageous for quarks to form pairs with non- 
vanishing total momentum 2q. This state, introduced for the first time in the sixties in the 
contest of electromagnetic superconductors, is known as Fulde-Ferrell-Larkin-Ovchinnikov 
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see 



2J for 



(FFLO) phase [21[. Its relevance in two flavor QCD has been discussed in 
a review). By virtue of the non- vanishing pair momentum, only a small region of the Fermi 
sphere is interested in the pairing phenomenon. This results in condensates that are smaller 
than the BCS ones. Moreover, the quark condensate in the FFLO phase is space-dependent. 
As a consequence, the translational and rotational symmetries are spontaneously broken 
and the spectrum of the low energy excitations is enriched by the presence of the respective 
Goldstone bosons, namely the phonons. As far as instability is concerned, the authors in 2^ 
have shown that, with two flavors, the instability of 2SC implies that the FFLO phase is 
energetically favored. However, the problem of the chromo-magnetic instability of the FFLO 
phase is under debate [3, • 

At densities relevant for the physics of the compact stellar objects the three flavors u, 
d and s can form FFLO pairs. In [2^ a first attempt to the study of the three flavor 
FFLO phase (at zero temperature) has been presented, based on a Ginzburg-Landau (GL) 
expansion of the pressure. The assumed pairing ansatz is 

3 

< AaC'y, i,fi, >= e"^' e,,-, (1) 

1=1 

with 

A7(r) = A/exp(2iqi-r) , (2) 

and the A/'s are the gap parameters. In the above equation 2qi represents the momentum 
of the Cooper pair. For values of the strange quark mass such that the FFLO phase is 
energetically favored with respect to the homogeneous phases, and among the different 
structures considered, it was found that Ai = and A2 = A3, q2 = qa is the configuration 
with the highest pressure. 

Recent interest has been devoted to the studvof the relevance of a lepton chemical 
potential on the phase diagram of Q, y l8ll |. This problem is directly connected 

to the phenomenon of the neutrino trapping, which could occur in the first cooling "era" 
of a proto-neutron star. In more detail, it is well known that the cooling of a neutron star 
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occurs via the emission of neutrinos 32, |33| (there is also a black body contribution due 



28| |. investigating the role of a non-zero 



to photon emission, but it is irrelevant for our scopes); if the temperatures are not lower 
than ~ 1 MeV then there is a spherical inner region in the star, the neutrinosphere, from 
where only a very small fraction of neutrinos escapes. The radius of the neutrinosphere 
(measured from the center of the star) depends on the temperature, and decreases as the 
star cools. Therefore the trapping interests regions closer and closer to the center of the 
star as its temperature lowers. As a consequence of the trapping there exists a (quasi)- 
conserved charge, the lepton number, so one can introduce a lepton chemical potential fii 
associated to it. If color superconductive quark matter is present in the cooling star, it 
would be interesting to investigate about the effect of the neutrino trapping on its structure. 
Recent works along this line show that a non-zero fiL has the effect to favor the 2SC phase, 
disfavoring the CFL phase js^ Isi | . 

In this short note we extend the results of 
on the three flavor FFLO phase of QCD. Since the neutrino trapping requires a non-zero 
temperature, we work at finite temperature in all this Letter (we should mention here that 
the critical temperature of the FFLO phase is expected to be lower than the homogeneous 
one Q)- 

II. THE MODEL 

We consider an electrical and color neutral system of massless u, d and massive s quarks, 
in /3 equilibrium with massless electrons and their neurinos. The lepton sector is described 
by the Dirac lagrangian 

£i = ilj{i^+ , (3) 

where we have collected the electron and the neutrino fields into the doublet = (e, z/), and 
the chemical potential matrix is fii = diag(/ie, fJ'u) = diag(— yUg + /i^, /i^); fiq is the chemical 
potential associated to the conserved electric charge of the system. From the lagrangian Q 
one derives the pressure js^ 

Pi = ^^3i dke\og{l + e'^) , (4) 



a=e,v 



where gi is the degeneration factor (that counts the spin degrees of freedom): it is equal to 
2 for electrons and 1 for neutrinos. One can check the correct normalization of p;, that in 



the limit T ^ becomes pi = fit/Un'^ + At^/247r2. 

Next we consider the quark sector. The quark lagrangian is 

= M + fi'yo)^ ~l G ^^^XaiJ ijYKi^ , (5) 

o 

where M°^^ = 6°'^ diag(0, 0, M^) is the current mass matrix; fi^^ is the matrix of the chemical 
potentials: they depend in general on /i (the baryon chemical potential), /ig, and fi^, /ig, 
related to the conserved color charges: 

where T3 and Tg are the usual SU{3) generators. Actually one should consider eight color 
chemical potentials, one for each conserved color charge =< ipTaip > js^; we have checked 
explicitly that only and ng can be non-zero. This motivates the choice in Eq. Q. 

The interaction term in Eq. (jSI) is a Nambu-Jona Lasinio inspired four fermion interaction, 
that mimics the one gluon exchange of QCD. Here G is a coupling constant, with dimension 
mass~'^; Xa are color matrices and a sum over flavors is understood. In the mean field 
approximation, after a Fierz rearrangement, the interaction term becomes 

-^6„^,e^^-^(C< A,(r) + c.c.) + {L ^ R) - ^Aj{r)A*j{r) , (7) 

where the assumed pairing ansatz is in Eq. ((T)). In getting Eq. ((Zj) we have neglected the 
chiral condensate: its effect is the dressing of the bare quark masses; we expect that at 
intermediate densities (namely fi ~ 500 MeV) the chiral condensate is small if compared 
to the quark-quark condensate. Therefore we can safely neglect the constituent u and d 
masses. As for the strange quark, we should write a gap equation for its constituent mass; 
it should to be solved simultaneously to the gap equations for the di-quark condensate. A 
similar analysis has been performed recently in ^ . For simplicity we do not solve the chiral 
gap equation and we assume the strange quark mass as an external parameter. 

We now discuss the other approximations used in the quark sector. First, we consider 
only the leading order effect of the strange quark mass, namely a shift in its chemical 
potential /i^ — * /i^ — M|/2/i. Second, to ensure color and electrical neutrality of the system, 
the chemical potentials related to the conserved charges have to satisfy the stationarity 
relations 

dp dp dp ^ 



where p = pi + Pquarks (see below, Eq. ffTUj) ). As in Ref. [28|, we work in the approximation 
= = 0. This assumption is justified because the phase transition from the supercon- 
ductive to the normal phase is second order. Therefore, near the phase transition we expect 
/i3, yUg not too different from their values in the normal phase, namely zero. In this way one 
has 

2 1 1 ^2 

fiu = fi+ ■^fJ'Q , fJ'd = fi- g/iQ , /is = /i - -/iQ - ^ • (9) 

Moreover we assume Ai = 0. This is justified because at T = and /i^ = one has 
fiQ = —M^/Afi. As can be inferred from Eq. Q this implies that the difference of the 
chemical potentials 6fids between d and s quarks is greater than Sfi^d and Sfius- As a 
consequence, the pairing between d and s is disfavored. This is true also if T 7^ and 
/ii 7^ (see below). Finally, the interaction contribution to the pressure is derived in the 
Ginzburg-Landau (GL) expansion. Keeping this in mind the pressure of quark matter is 



Pq = Po- 



Al + ^At + ^Al + ^At + ^AlAl + 0(A«) 



(10) 



2 ^ 4 z 2 ^ 4 ^ 2 
where po is the pressure of the normal phase. In Eq. pO|) the coefficients are given in terms 



of the functions 



a{q5ix,T) 
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as follows 



' GAtt^T^ Ait ^ \2 2-kT 



(12) 



a2 = a{q2, 03 = "(gs, 7: , -l ) 



P2 = P{q2 



2 

/ili /is 



2 

/id — f^u 



T) 



moreover 
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S/i^ 1 



dyidy2S{l - yi - y2) / rfzRe^/'^^M - 



(13) 



In the above relations fi is the baryon chemical potential as introduced in Eq. the 
functions tp^^-^z) are defined as ^^''-derivatives of the Euler ip{z) function, where ip{z) = 
T'{z)/T{z). In Eq. (fT^ we have introduced the function 

< /iu /is /id /i?i ,1 , I 

^ = 1/1 n ?/2 ^ h 2|l/iq2 + Z/2q3| • 



The angular integrals in Eqs. (fTTj). (fT^ and (fT^ can be done exactly, but their expressions 
are uninformative so we prefer to leave them in implicit form. As for the integrals in yi, 
1/2 in Eq, (fT^ . one is performed by the 5(1 — j/i — 7/2) and the remaining integral can be 
performed numerically. 

dp 

The magnitudes of the wave vectors q/ are determined by the variational condition — — = 

Qoqi 
, in this paper we use this relation at the lowest order, — — = 0. If T = this 

oqi 

approximation leads, as in the two flavor case, to the well-known relation qj = 1.1997 



III. RESULTS AND CONCLUSIONS 



In this section we discuss our results; the coupling constant G can be eliminated in favor 
of Aq, the gap parameter of the homogeneous three flavor superconductor at T = and 

= Q- We show the results obtained for Aq = 25 MeV, fi = 500 MeV and T = 0.1 Aq 
(for different values of the parameters we get qualitatively similar results). In Fig. C] we 
show the electrical chemical potential fig that satisfies the stationarity condition (jH)) as a 
function of M^/fi for four values of the lepton chemical potential fiL- We notice that for a 
fixed value of the strange quark mass, increasing fi^ results in the decreasing of This 
has important consequences on the pairing. Indeed one has from Eq. Q 

For fiL = one has fig = —M^/Afi; therefore the mismatch between the u and the d Fermi 
surfaces is the same of the u and the s one. For fiL ^ this is no longer true: in particular 
from Fig. ^one gets fJ^u — fJ's > f^d — fJ'u- This favors the pairing in the u — d channel while 
disfavors the u — s pairing. 

This is better expressed in Fig. |2l where we show the results for the gap parameters A2, 
A3 as a function of M^//i, for /i^ = (left panel) and /iz, = 200 MeV (right panel). In 
the latter case we observe that for low values of the strange quark mass the gaps A2 and 
A3 have comparable magnitude. This means that u ~ d and u — s pairing are both active 
and the FFLO state is effectively a three flavor color superconductor. Increasing Ms we find 
the window 100 < M^/fi < 120 MeV where both u — d and u — s pairing are active, but 
A2 < A3 meaning that the latter is disfavored if compared to the former. At / /i ~ 120 
MeV a second order phase transition takes place to the state with A2 = and A3 7^ 0, that 
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FIG. 1: Chemical potential /ig as a function of M"^ / ^ at T = O.IAq. The cases = 0, 100, 200 and 300 
MeV are represented respectively by the solid, long-dashed, short-dashed and dot dashed lines. 




FIG. 2: Left panel: A2/A0 = A3/A0 against M^/^, evaluated for = 0. Right panel: A2/A0 (dashed 
line) and A3/A0 (solid line) as a function of M"^ / ji, computed for = 200 MeV. In both pictures T = O.IAq. 

is a two flavor FFLO state with pairing in the u — d channel This is in perfect agreement 
to what is found in the homogeneous case 3^ 0|. Finally, for M^//i ~ 185 MeV there is 
a second order phase transition to the non superconductive phase. We find also that the 
larger is the value of /ij^, the larger is window of M^/yU where A2 = and A3 7^ 0. 

In conclusion, we have investigated the role of a lepton chemical potential on the structure 
of the three flavor FFLO phase of QCD. We find that a non- zero strongly favors two 
flavor pairing, in the channel u — d, while disfavoring the pairing u — s. This is a result of 
imposing electric neutrality in the system. 

This study could be relevant for the cooling of a proto-neutron star. If one suppose that 
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normal (that is, non superconductive) quark matter is present in the star, then the neutrino 
emissivity will be dominated by direct URCA processes S^Q; denoting by e^^'-^^{T) the 
neutrino emissivity of normal quark matter in absence of trapping at the temperature T, 
the effect of trapping is an exponential suppression of the emissivity itself jsTj 



e.(r, T) = exp [-j^) x e,^«^^(T) ; (14) 

here ey{r,6,T) denotes the emissivity with the effect of the trapping included; l{r,6) is the 
distance from the creation point of the neutrino to the surface of the star (more precisely, its 
projection along the 2;-axes). Finally, A(T) is the neutrino mean free path at the temperature 
T. The exponential factor in Eq. (fT^ takes into account the probability that a neutrino 
created at a distance r from the center of the star can leave the star in the direction defined 
by the angle 9. One computes the total luminosity for neutrino emission from the star by 
averaging Eq. (|14j) over all neutrino directions 9 and integrating over all distances r up to 
the star radius. 

We turn to color superconductive quark matter. Un-trapping occurs for temperatures of 
order of 1 MeV (for higher temperatures, neutrinos are trapped); this imply that trapping 
effects should not be neglected in emissivity computations, unless one considers the final 
cooling epoch of a neutron star. Investigations on the effects of color superconductive quark 
matter on the cooling evolution of a neutron star, neglecting neutrino trapping, have been 
performed in 



39| (see also references therein). It will be interesting to see the effect of 



the FFLO phase on neutrino emissivity, if such a phase is actually present in the core of 
compact stars. 
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